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Abstract 
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exponential generating function of correlators is calculated explicitely. The calcula- 
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terminant formulas for the norms of Bethe eigenfunctions due to Richardson and 
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1. Introduction 

The most challenging problem of the theory of quantum integrable systems, 
after describing the spectrum and eigenfunctions, is finding a reasonably efficient 
expression for the correlators. Speaking about the systems lying in the realm of 
Bethe ansatz (R-matrix, quantum inverse scattering) method, two complementary 
approaches are presently known. One of them is based on studying analyticity and 
bootstrap equations for factorized S'-matrices and formfactors of infinite-volume 
quantum field models and is useful in the study of short-distance behaviour of cor- 



relators (see the review papers fl4] , g). Another one, summarized in the book ||, 
is, on the contrary, efficient in describing the long-distance asymptotics and is based 
on meticulous study of the structure of Bethe eigenfunctions in the finite volume. 

Whereas the mathematical structures involved with the former approach are 
now recognized quite satisfactorily (quantum Knizhnik-Zamolodchikov equation, 
vertex operators, representations of quantum groups) which has resulted in the rapid 
progress in this subject, this is not the case for the latter approach. After tour de 
force of 1982 by V. Korepin JF], who had managed to overcome the enormous combi- 
natorial difficulties in handling the cumbersome expressions for correlators produced 
via algebraic Bethe ansatz, very little progress has been achieved in mathematical 
understanding of his technique. One should mention || where Drinfeld's twists were 



applied to simplify Korepin's calculations, and ]10| where an alternative derivation 



of the determinant formula for the Bethe eigenfunctions based on the semiclassical 
asymptotics of the solutions to the Knizhnik-Zamolodchikov equation was proposed. 

The present paper has grown from an attempt to understand and to simplify 
Korepin's derivations using a simple toy example. As such, I have chosen the sl 2 
Gaudin model [T], |3| (a degenerated case of integrable XXX spin chain). The choice 
was motivated by the fact that Gaudin model had already played this role once. 
The determinant formula for the Bethe eigenfunctions of Gaudin model obtained 



by Gaudin, who, in turn, had benefited from an earlier paper by Richardson |TT 
allowed to Gaudin, McCoy and Wu to conjecture in [Q the similar formula for the 
XXX spin chain which was finally proved by Korepin [^]. 

The intimidating combinatorial complexity of the algebraic expressions for cor- 
relators resulting from algebraic Bethe ansatz makes it natural to take advantage of 
the common combinatorial wisdom which says that using generating functions can 
help to reduce the complexity drastically. As shown in the present paper, this simple 
idea works quite effectively in case of the Gaudin model. A simple explicit expres- 
sion for the exponential generating function of correlators is obtained from the Gauss 
decomposition for SL(2) loop algebra. From the generating function the formulas 
for particular correlators are derived quite easily, including the Richardson-Gaudin 
determinant formula for the norm of Bethe eigenf unction. 

The paper is organised as follows. In Section 2 we introduce the sl 2 Gaudin model 
and present the known results for its eigenfunctions, in particular, the Richardson- 
Gaudin fomula. In Section 3, we describe a new formula, called A-representation 
for correlators which we prove in the next, 4th section using the generating func- 
tion technique. The idea of the proof is explained first on the elementary finite- 
dimensional example of sl 2 Lie algebra. The A-representation is used in Section 5 to 
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derive the Richardson-Gaudin determinant formula. In the concluding, 6th section 
we discuss the obtained results and the prospects of their generalisation. 

This work has been started during my stay at the Department of Mathematical 
Sciences, the University of Tokyo, when I was preparing a course of lectures for 
graduate students on Gaudin model and was continued in Research Institute for 
Mathematical Studies, Kyoto University. I am grateful to UoT. and RIMS for 
support and hospitality. 



2. Description of the model 

The aim of this section is mainly to fix the notation. The results given below 
are taken from from |j[ |13| . 

Let e, f, f) be the generators of the Lie algebra s^: 

[e,fl = f) [M = 2e [M = "2f, (2.1) 

and let a highest weight module labelled with the complex number (highest weight) 
A be generated by the vacuum (highest vector) |0) and the relations 

e|0)=0, f>|0) = A|0). (2.2) 

We shall need also the dual module generated by the dual vacuum (0| and rela- 
tions 

(0|f = 0, <0| f) = A <0| , (0|0) = 1. (2.3) 
The corresponding value of the Casimir operator 

2K = ef + fe + i() 2 (2.4) 

is then 2K = A(A + 2)/2. 

The commuting Hamiltonians Sj of the Gaudin model are defined in the tensor 
product of D h.w. modules (we label the corresponding generators with the index j 
and put |0) := |0) 1 <E> ■ ■ • <S> |0) D ) and are given by the expression 

s . = E ^±i^±lMi, J = 1,...D (2.5) 

where Zj are complex constants. It is more convenient, however, to express Hj in 
terms of a generating function. To this end, consider the one-parameter operator 
families 

= F(u) = f:-^-, H{u) (2.6) 

j=l " ~? j=l " ^3 j=l " ^3 

and define the generating function t(u) as the result of replacing c, f, () in the Casimir 
( |2.4j ) respectively with E(u), F{u) and H(u): 

2t(u) = E{u)F{u) + F{u)E{u) + \h 2 (u) (2.7) 
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Using (|2.6p and Q2.5|) one observes that 



The fundamental property of the operators ( |2.6| ) is that they form the highest- 
weight module over the infinite-dimensional Lie algebra f||, [13| 

[£(«), £(»)] = [F(u), F(t>)] = [//(«), //(b)] = 0, (2.9a) 
[g( B ),F(.)l = - g( " ) - g( ' ) . (2.9b) 

[g M ,F W ] = -2 g W- g W , (2.9c) 

u — V 

mu),F(v)]=2 F M-_y , (2.9d) 

characterised by the vacuum |0), 

E(u) \0) = 0, H(u) |0) = A(w) |0) , (2.10) 

dual vacuum (0| 

(0|F(u)=0, (0|^(u) = A(«)(0|, <0 | 0> = 1, (2.11) 
and the scalar function (highest weight) X(u) 

A(tt) = f;-^-. (2.12) 

3=1 " ~7 

A direct consequence of the relations ( |2.9|) is, in particular, the commutativity 
of t(u) 

[t(u),t(v)] = (2.13) 

from which the commutativity of the hamiltonians Hj follows immediately. 

The commutation relations Q2.9|) can be also written down compactly in the 
so-called r-matrix form |3|, [12|, [L3| which we, however, shall not use here. 

The theory of Gaudin model is expected to give solution to two main problems. 
The first one is to determine the eigenvalues and eigenvectors of the commuting 
Hamiltonians Hj ( |2.5| ). 

The solution is given in terms of Bethe vectors which are defined for any finite 
set V of complex numbers as 

|V> := J] F(v) |0), VCC, |V|<oo. (2.14) 
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Theorem 1 [73| / The vector |V) is a joint eigenvector of commuting operators 

Ej or, equivalently , of t{u) if and only if the parameters v G V satisfy the Bethe 
equations 

X ( v ) = E ——n VueV - ( 2 - 15 ) 



, , v — v 

The corresponding eigenvalue t(u) of t(u) is then 

T(u) = ~\ 2 (u)-~d u \(u) (2.16) 

where 

\(u) := X(u) - E (2-17) 



The spectrum of the model having been determined, the next fundamental prob- 
lem is to calculate the correlators. The correlators which we are going to study are 
labelled by three finite sets U, V, W C C and are defined as 

C(U,W,V) := (U\ I] H(w)\V) (2.18) 
wew 

where {U\ is the dual Bethe vector 

(U\ := (0| J] E(u). (2.19) 

In principle, one could calculate C(U,W,V) moving the operators H{w) and 
E(u) to the right through F(v ) with the use of the formulas 

F(u)|V) = |VU{u}), (2.20a) 
H(w) |V) = (\(w) - J2 -M IV) + E 77^- IV \ W U {w}) , (2.20b) 

E{u) |V) = £ --i- (\{u) - X(v) - £ 2 f-^-j - ) |V \ {v}} 



+ E V( mIV\K^'}uW) (2.20c) 

which follow from the relations ( |2.9| ) and fl2.10| ). For example 

C({u},$,{v}) = - X{u) ~ X{V \ u + v. (2.21) 

u — V 

This strategy was adopted in where the formulas similar to (|2.20|) were used in 
a more complicated case of XXX spin chain. Such straightforward approach leads, 
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however, to extremely tedious calculations, and it is the aim of the present paper to 
simplify the derivation using the method of generating functions. Still, one simple 
observation can be made right now: the correlator C(U,W,V) is a polynomial in 
A(it), A(t>), \(w) with the coeffitients rational in u G U, v G V, w G W. 

It is important to stress, following |7) , that for the derivation of the formulas for 
C(U,W,V) the concrete form ( |2.12j ) of the function A(w) is not important, it can 



be considered as a formal parameter. All what matters is analyticity of X(u) at the 
points of U U V U W and the relations (Q, ( P^UQ and (gTg ). 



We conclude this list of known results with an important determinant formula 
for the norm of the Bethe eigenfunction. 



Theorem 2 Let V = {v\, . . . ,vn} and Vj satisfy the Bethe equations ( 2.15] ). Then 

(V | V) = (-l fdetM (2.22) 
where the matrix Ai is given by 

Ma = d Vj X( Vj ) + £ _ 2 (2.23a) 

k^j \ u i °k) 

M 3k = — — -, j^k (2.23b) 



The formula ( 2722 ) was derived first by Richardson |TlJ for a model which can 



now be recognized as a very degenerate case of Gaudin model. In || Gaudin men- 
tioned that the norms for Gaudin model are given by the same formula but had not 
explained how to generalise Richardson's derivation. To make up the deficiency we 
give a complete proof of ( [2.22j ) in section 5. 

3. A-representation for correlators 

In this section we give a new explicit formula, which we call A-representation, 
for the correlator C(U, W, V). Throughout this paper we suppose, unless otherwise 
stated, that the sets U, V and W do not pairwise intersect. 



From fl2.20|) it follows immediately that C(W,W, V) is a polynomial in A(x), 
x G IA U V U W with the coefficients rational in x. However, the structure of 
C(U, W, V) can be simplified after introducing a proper notation. 

For any finite set X C C define the function 

A(X) := V — ^ r . (3.1) 



For example, 



A({x}) = A(x), A({x 1 ,x 2 }) = A(Xl) A(X2) . (3.2) 

X\ — %2 
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Note that, supposing X(u) to be analytic inside a contour T circumscribing coun- 
terclockwise the set X, one can rewrite (|3.1| ) as the contour integral 



AW = 1 / ^ )dZ r , *cC + . (3.3) 

xex 

Our aim is to express C(U, W, V) as a polynomial in A. 
Let |-| be the cardinality of a set, and suppose 

\U\ = |V| = N, |W| = M (3.4) 

(obviously, C(U, W, V)=0 if \U\ ^ |V|). Let us define a coordinated partition V of 
the sets W, V and W as a set of triplets P = (U P , V P , Wp), U D U P ^ 0, V D V P ^ 0, 
W D Wp (note that Wp, in contrast with U P and Vp, is allowed to be empty) such 
that 

VP, P' EV : |Wp| = |V P | > 0, I Wp| > 0, 

p^p' UpDUp, = 0, Vp n v P , = 0, Wp n w P > = 0, 

U Wp = u, \jv P = v } (J Wp c w, 

p e p p e p p e p 

E |Wp| = |w| = E \Vp\ = M = n, E \w P \ < |w| = m. 

Pep Pev Pev 

Theorem 3 The expression for the correlator C(U, W, V) is given by the formula 
C(U,W,V) 

= (-1)* E f II W (\Up\ ~ !) ! (2 |Wp|) |w ^ A(Wp U Vp U Wp)) 

p \pg-p / 



x 



n a h I ( 3 - 5 ) 

\wew-p 



where 

Wp := W \ (J Wp. (3.6) 



PeP 



For example, 

C({u}, {w}, W) = -2A({u, «}) - A({u, v})A({w}). 

It is remarkable that C(JA, W, V) is a polynomial in A with integer coefficients, 
that is all the rationality is hidden in the definition ( |3.1| ) of A. 

To describe the norms of Bethe eigenfunctions we shall need a particular case of 
Q when W = 0: 

(U\ V) = (-l) iV E II |Wp|!(|Wp|-1)!A(W p UVp) (3.7) 
v Pev 
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where P runs over the sets of pairs P = (Up, Vp), U D Up ^ 0, V D Vp 7^ 0, such 
that 

VPEV : \U P \ = \V P \ , 

P,P' eP, P + P' =^ U P n Up, = 0, V P n v P > = 0, 

U Up = u, (J v P = v, 

PeP PeP 

X; |W P | = \U\ = X |Vp| = |V| = N. 
PeP PeV 

For example, 

(0 | 0> = 1, (u\v) = -A({u,v}), (3.8a) 
(uiu 2 I ui-ua) = A({wi, wi})A({u 2 , v 2 }) + A({«!, w 2 })A({m 2 , ui}) 

+2A({ui,U2,«i,u 2 }). (3.8b) 



Using the integral formula ( |3.3| ) for A(X) it is easy to extend by analyticity the 
formula (|3.5|) for C(U,W,V) to the case when the sets U, V, W intersect or even 
have multiple points (being thus not sets but divisors). We leave this exercise to the 
reader. 



4. Proof of the Theorem [3| 

As mentioned in the Introduction, we are going to circumvent the combinatorial 
difficulties of the direct approach based on the identities (|2.20|) by calculating a 
generating function of correlators C(U, W, V). 

Let us illustrate the idea on an elementary finite-dimensional example of the s/ 2 
Lie algebra Q. The analog of C(U, W, V) is the expression c(k,j) := (0| e^f |0) 
which can easily be estimated as 

c(k,j) = (-k) k (-X) k (X - 2ky = (-l) k k\ (-A) fc (A - 2k) j (4.1) 

where the Pochhammer symbol is defined as 

(a) k := a(a + 1) . . . (a + k - 1). (4.2) 

On the other hand, let us introduce the group elements 

S e :=e £t , W„:=e"&, T v := e^, (4.3) 

parametrized by coordinates e, 77, (p and consider the exponential generating function 

oo oo k j k 

<o| 4^ |o> = E E £ Jr|- c (M. ( 4 - 4 ) 

Notice now that if we manage to change the order of factors 

S £ Ti. ri J-' !fi = J r !fi iTi, ri iS e i , (4-5) 
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then the generating function is easily calculated: 

(0| ^HrfEe' |0) = (0| Hrf |0> = e"' A (4.6) 

To establish the identity (|4.5|) it is sufficient to consider the fundamental (2- 
dimensional) representation of Lie algebra sl 2 

-(Si), •-(ii). »-(:$)■ <-> 

and the corresponding Lie group SL 2 

The identity (O) expresses then equivalence of left and right Gauss decomposi- 
tions in SL2. Multiplying 2x2 matrices we find from (|4.5| ) the values of the primed 
parameters: 



= + £ ' = vT " -v V = (4-9) 

e v + eipe e^ + eipe 



whence 

(0| E S H^F V |0> = (e" + eye^f = e Ar? (l + £^r 2 ")\ (4.10) 
Applying the binomial expansion formula 

(i _*)- = £; Mi z * 

k=o K - 

to ( |4. 10|) one recovers the result (|4.1| ). 

The same Gauss decomposition trick works also in the case of the infinite- 
dimensional Lie algebra (|2.9| ). 

Let C + be a simply connected bounded open domain on the complex plain hav- 
ing smooth boundary T which is homeomorphic to a circle S 1 . Let C_ be the 
complement domain: C_ = C \ (C+ U T). Suppose that the contour T is oriented 
counterclockwise. Suppose also that U, V, W C C + and that X(z) is analytical in 

c+. 

Let e(x), ip(x), r](x) be some smooth functions on T. Define the generators of 
the Lie algebra A as 

E £ = ^- [ e(x)E(x), H v = f v (x)H(x), F v = / <p(x)F(x), 
2m Jr 2m Jr 2m Jr 

(4.11) 

To describe the commutation relations between them we need to introduce the 
decomposition 

if)(x) = i/j+(x) + if)-(x) 
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of any function ip(x) on the contour T into the parts ip±(z) analytical respectively 
in z G C±. It is supposed that ip-(z) — ► as z — > oo. The projections are given by 
the Cauchy integral 

. . 1 r ib(x)dx „ _ 

^ ± z :=±— / , xGl\ 2GC ± . 

Note the formula for the Hilbert transform defined as the singular integral 

2m Jr x — y 2 

regularized in the sense of the principal value. 

Now we can write down the commutation relations for the generators (|4.11| ) as 

[E E (i),E £ (2)] = [H v (i) , H v (2)) = [F v (i), -Fp(2)] = 0, (4.12a) 

[E e ,F v \ = H e+v+ - e _ v _, (4.12b) 

[H v ,E e ]=2E v+s+ _ ri _ E _, (4.12c) 

\Hni Ftp] 2F r j +! p + _ r i_ t p_ , (4.12d) 



The above formulas suggest the decomposition, see p2| , of A into two mutually 
commutative subalgebras A = A + + A- generated, respectively, by E £± , H v± , F ip± : 

E e = E £+ + E e _ , iJ^ = H v+ + H v _ , F v = F v+ + F v _ . 

The commutation relations for the subalgebras A± are those for the s/2 loop Lie 
algebra 

[E £ m , Eg] = [# « , # g)] = [F^j , F g)] = 0, (4.13a) 

[F £± , F^J = ±H e±ip± , [H V± ,E £± ] = ±2E V±£± , [H v± , F v± \ = ^2F v±(p± , (4.13b) 
and allow the fundamental representation as 2 x 2 matrix-valued functions of x G T 

E £± = ±£ ± (x) (° iy H n± = ±T] ± (x) ^ J _^ ^ , F v± = ±(p ± (x) ( J q ) • 

with the pointwise commutator. 

Now everything is ready to reproduce the same Gauss decomposition trick which 
we saw working in the finite-dimensional case. 

Introduce the group elements 

£ e _j_ '— exp E e j s _ , H"rj± - = sxp-ffjy., J~ p± '.— exp-Fj^, 
having the fundamental representation 

/ 1 ±e ± (x) \ ( e±^) \ ( 1 \ 

Ce± ~ ^ 1 j> n v±- y e^±w J ' ^-y ±(p ± ( x ) 1 y • 
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The identity 

£e±H r}± J : l p ± = ^>' ± W^£ £ ' ± (4-14) 

is established by the same calculation as in the finite-dimensional case 

e ±V± = e ±r,± + £±(f±e TV± ; (4. 15 ) 

± e ±? ?± + e±ip±e^ r i± ' V± e ±? ?± + £± ^ ±e ^± ' 1 ' ' 

one should rememeber only that now e±, etc. depend on parameter x G V. 

The formula for the generating function of correlators follows then immediately: 



= (Fip' + T~L V ' + £e' + ) (•?>'_ 7~Lr{_ £e'_ 

= (H v > + ) (7V) = exp ^ | r (4.17) 

where 

rf = rf + + rf_, rj' ± = ±ln(e ±v± + e TV± e±tp±). (4.18) 
Now we can prove Theorem |[ Put 

e{z) := £ -r-, ^) : = E *>(*) : = E (4-19) 

where e u ,r] w ,ip v G C. Then, since by assumption U C C + , 

„ 1 /■ v-^ e u E(z)dz ^ . , 

£i '=« p 5s/5-THr = 6,p 5 8 * B w- (420) 



and similarly for i^, i 7 ^. The correlator C(W, W, V) is given then by the coefficient 
at 



\uau J Wv / \weW / 



in the expansion of the generating function (|4.17|) in powers of e u , r) w , <p v . 
To perform the expansion, note first that 

e+(z) = 0, e.{z) = e{z) (4.21) 

and similarly for <p(z), T}(z). Respectively, rf + = and 

!]=!]_ = - ln(e~ v + epe v ) =r)- ln(l + eipe 2v ) 

= ^-EE ( i ( 7 W, (4-22) 

k=lj=0 K 3- 

where we omitted the argument z of e(z), f](z), <p(z). As it was explained above, we 
are interested only in the terms linear in e u , tp v , r] w . For example, the contribution 
of e k is 

*l E (4-23) 

Wcw new ^ " 
|W'|=fc 
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Respectively, the relevant terms in rj'(z) are 



Err+ E E E (l^'l - !) ! ( 2 l^'l) |w '' 



ygW ^ W u'cu v'cv W'CW 



|w'| = |v'| v'#8 



X 



J J g tt J J (fv J J <]w 

ueU' z ~ u „eV Z ~ V lugW z ~ w 



(4.24) 



Integrating the result with A(z) and using the equality ( |3.3| ) we obtain finally for 
the relevant contribution to (1/2-ki) J X(z)i]'(z)dz: 

E vwKw) + E E E |w / |!(|w / |-i)!(2|W|)i w 'i 



w'cw v'cv W'CW 
|W|=|v'| vV0 



X 



n e„ ) ( n v« ) ( n ^ ) a ( w ' 

K uew ) \vev ) \wew J 



U V U W) (4.25) 



Exponentiating ( f4.25|) and retaining the term containing the factor Y[£ u <^ v r] w we 
finally obtain (3.5). 



5. Richardson's determinant formula 

For the rest of this paper we consider only the case W = 0. To prove the Theorem 
H, we derive first a recurrence relation for the correlators C(U, 0, V) = (U | V). 
Note the following property of the function A 

A( * uy) -s£ ji(«- A srp!,-^) <5 - i) 

which holds for X fl y = and can be deduced easily from the definition ( p.l| ) of A. 

Using (|5.1|) one can rewrite the A-representation Q3.7|) for (U | V) as a polynomial 
in the two-point correlators A({u, v}) with the coefficients rational in u, v. Let us 
denote the polynomial in question as C(U, V; £) considering £(u, v) = A({w, v}) as a 
functional parameter: 

(U\V)=C(U,V;A({u,v})). 
Richardson has found an alternative description of the function C(U, V;£). 

Theorem 4 The function C(U, V; £) possesses the following properties: 

1. C(U, V; £) is a polynomial in £(u, v), u G U, v 6 V. 

2. Each of the £(u, v) enters C(U, V; £) linearly. 

3. The polynomial C(U, V; £) has no free term, that is 

C(W,V;0) = 0. 
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4- Recurrence relation. Let U = U U {u} and V = VU {v}. Then the coefficient 
at £(u, v) in C(U, V; £) is equal to —C(U, V; £) where 

I(u, v) := £(u, v) + ^ =v ( 5 - 2 ) 

(u — u)[v — V) 

5. Initial condition 

C({u},{v};£) = £(u,v). 



Obviously, the properties 1-5, if true, characterize the function C(U, V; £) uniquely. 

Richardson himself has proved the above theorem by a straightforward calcu- 
lation for a very degenerate case of Gaudin model. Our proof, valid for general 
Gaudin model, is based on the A-representation for (U | V) and is much simpler due 
to using the machinery of generating functions. 

The statements 1, 2 and 3 of the Theorem follow directly from the A-representation 
( |3.7p and ( |5 . 1| ) , the statement 5 follows from ( |2.21| ) and (|3.2| ). The remaining state- 



ment 4 (recurrence relation), in principle, can also be derived directly from the 
A-representation which leads, however, to lengthy combinatorial calculations. As 
in case of the theorem ||], using generating functions again allows to simplify the 
derivation, . 

Putting W = in fl4.25|) one obtains the relevant (that is linear in e u , (p v ) 
contribution to (l/27ri) / X(z)rj'(z)dz: 

e e (-i) 1 "' 1 m - 1)! ( n £u) [ n <*) hw u vy ^ 

u'cu v'cv \ueW J \veV J 

|w'|=|v'| v'#0 



Let us find, first, how the substitution £ i— > £ (|5.2j ) can be expressed in terms of 
the expansion (|5.3j). From ( p.l| ) it follows that replacing all A({u, v}) with ( |5.2j ) is 



equivalent to replacing all A(W U V) in (|5.3|) with 




Substituting ( p.4| ) into ( |5.3| ) we conclude, similarly to the derivation of ( |4.24| ) 



that the terms added to (l/2iri) J \(z)rj'(z)dz coincide exactly with the expansion 
of 

-21n(l +e(u)(p(v)) 



where e(z) and (p(z) are given by ( |4.19| ). 



Exponentiating, one obtains that the substitution £ i— > £ is equivalent to multi- 
plying the generating function {E e T^) by the factor 



exp {-2 In (1 + e(u)(p(v))} = (1 + e{u)^{v)y 2 . (5.5) 
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Using the binomial series 

oo 

(i+r 2 = E(-i) fc (^+i)^ 



k=0 



one can expand the expression ( |5.5|) in powers of e u , tp v (see again the derivation of 
(4.24) from Q4.22p ) and retain the relevant terms obtaining finally 



E E(-i) |w ' l+1 OT + i) ! M! ( II ttztA (n^i- M 

U'cU 

\u<\=\v>\ 



u'cm vcv \u&A' u u V«eV v v 



Now we return to the expression ( |5.3|) for the relevant part of (l/2m) J X(z)r)'(z)dz 
and replace the sets hi by hi and V by V. In the resulting expansion 

e e (-i) 1 "' 1 w (iwi ( n f n a ( w/ u v) (5.7) 

wcu vcv \ueW J \veV J 

|w'|=|v'| vV« 

four types of terms can be distinguished. 

1. The terms containing no u, v reproduce obviously the unperturbed expansion 
Q. 

2. The terms containing u but not v. 

3. The terms containing v but not u. 

4. The terms containing u &W — hi" U {u} and v G V — V" U {w} 



w e e (-ir i+i aw"i+i)! iwi! n e« n ^ a^wur^) 

m"cw v"cv VugW" / VtieV" / 

| W »|=|V"| 

(5.8) 

(note that here hi" and V" are allowed to be empty. 

Obviously, only the terms of type 4 produce A({u,v}) via the expansion (|5.1|). 
Noticing that the coefficient at A({u,v}) in the expansion (|5.1|) of A(W'UV' U{u,v}) 
equals 



iK^r 1 n( 



V — V ) 1 



and exponentiating ( |5.8|) we obtain that the coefficient at e~<p^A({u, v}) is given 
exactly by the expression ( |5.6|) which proves the theorem |j. 

The recurrence relation allows to give for C(hi, V; €) a completely different repre- 
sentation. 

Theorem 5 (Richardson) Let us fix an ordering of the sets hi = {m, . . . ,Un} o,nd 
V = {v\, . . . , vn}- The function C(hi, V; £) can be expressed as the sum of N\ deter- 
minants 

{-1) N C{U,V;£)= E detAT (5.9) 
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where the sum is taken over all permutations a of indices {1, . . . , N}, and the ma- 
trices M. a are defined as 



M a - -- 



[Uj, v aj ) + 2 i u j ~ u f) 1 ( v o 

j'¥>3 



-2{uj - u k ) 1 (v CTj - v ak ) l , j ^ k 



(5.10) 
(5.11) 



Proof. It is easy to verify (see |TT|) that the right-hand-side of ( p.9| ) satisfies 
the same conditions 1-4 of the Theorem 1 which determine uniquely C(U, V; £). The 
theorem |2] is obtained now as a simple corollary. Indeed, if U = V and the parameters 
Vk satisfy Bethe equations ( 2.15Q then in the sum (|5.9|) only the term corresponding 
to the identical permutation is non-zero, and the formula for the norm of Bethe 
function simplifies to ( |2.22| ). For the details see again |]TT| . 

For a completely different proof of the formula (|2.22| ) based on the semiclassical 
expansion of the solutions to the Knizhnik-Zamolodchikov equation, see [1C 



6. Discussion 

On a simple example of the Gaudin model we have demonstrated that the calcu- 
lation of polynomial correlators simplifies drastically when one makes use of expo- 
nential generating functions. An intriguing and so far open question is if a similar 
simplification is possible in case of XXX or XXZ magnetic chains, when the Lie 
algebra (|2.9|) is replaced by the yangian 3^[s/ 2 ] or, respectively, quantum Lie algebra 
U q [sl-2\. It is natural to expect that the complicated formulas for polynomial cor- 
relators due to Korepin || |7j can be somehow simplified by appropriate choice of 
generating functions. It would be also interesting to find in this case an analog of 
our A-representation and to find an interpretation of Korepin's 'dual fields' in terms 
of generating functions. 

As a simple introductory exercise, one can consider the finite-dimensional quan- 
tum Lie algebra U q [sl<2\ and reproduce the derivation given in the beginning of the 
section 4. The necessary formulas concerning the Gauss decomposition for U q [sl2\ 
can be found in 0. The main new feature in comparison to the Lie algebraic case 
is that the parameters of expansion e, rj and <p become noncommuting. 
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